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Abstract 



We study the non-abelian, gauge invariant, screening potential between two charges 
in a quark-gluon plasma. Using Braaten and Pisarski resummed perturbation theory 
in temporal axial gauge, we find a repulsive power screened potential which behaves as 
~ 1/r 6 at large distance. This questions either the physical meaning of the Debye mass 
or the known pragmatic treatment of the temporal axial gauge propagator. 
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1 Introduction 



Charges in a medium are understood to be screened. At large distance, this screening 
is usually taken to be an exponential [|lj] with the screening length given by the inverse 
of the Debye mass m. In QED, according to linear response theory, this screening is 
related to the two-point function of the photon. This is due to the field strength being 
linear in the vector potential. The small momentum behaviour of this function gives 
rise to imaginary poles which leads to the exponentially screened potential, from which 
the Debye mass derives its physical significance of inverse screening length. 

In the last few years, the plasmon puzzle or the gauge dependence of the gluon damp- 
ing constant has made it clear that it is necessary to employ resummed perturbation 
theory of Braaten and Pisarski || whenever momenta much less than the temperature 
T of the system under study are involved. Since the long distance behaviour of the 
screening potential and hence the two-point function is required, it is necessary to use 
this resummed theory. 

Unlike QED, in QCD only the leading term of the static time-time component of 
the gluon self energy, i.e. the Debye mass squared, calculated in the framework of this 
resummed theory, is gauge independent 0. Therefore it is not clear that the Debye mass 
is a physical quantity in the non-abelian theory. By working out the leading correction to 
the Debye mass in the resummed theory, Rebhan || f|, [| argued that in covariant gauge, 
it was possible to have a gauge parameter independent and hence gauge independent 
Debye mass by defining 

m 2 = lim n 00 (0,q) , (1) 

instead of 

m 2 = lim n 00 (0,q) . (2) 

That is one obtains m self-consistently from the pole of the propagator. 

The leading correction to this gauge independent Debye mass was found to be 

5m 2 ~ gm 2 ln(l/g) . (3) 

The ln(l/g) is due to the introduction of the magnetic mass as a infrared cut off in the 
transverse propagator. As a result, the correction is rather large. The new definition Eq. 
(0) is possible because the gauge parameter dependent term in the leading correction in 
covariant gauge is proportional to q 2 + m 2 . However, it is not clear why m found this way 
in a gauge such as covariant gauge can be readily labelled as the Debye mass: the latter 
is generally understood to be the inverse screening length in a Yukawa type screening 
potential between charges and in covariant gauge, the potential is related to more than 

2 In the bare theory B, in addition to the gauge invariant leading term, there is a linear term in the 
modulus of the external momentum which is also gauge independent. 
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just the longitudinal propagator. On the other hand, the temporal axial gauge (TAG) 
is the only gauge in which we have this simplification. Baier and Kalashnikov have 
recently shown, by some unique construction, that the gauge-invariant definition of the 
Debye mass might not extend to TAG: instead of being exponentially screened, the 
potential is repulsive with power screening ~ 1/r 6 in TAG 0. It then seems appropriate 
to investigate this problem in TAG in the framework of Braaten and Pisarski resummed 
theory. 

The choice of TAG at finite temperature may be controversial, because it is known 
to have problems with the double pole when po is zero. Nevertheless, TAG has the 
advantage that one can relate unambiguously the potential between two charges to the 
two-point function. Also it has been shown in non-resummed theory || that the same 
result for n 00 (0, q — > 0) could be obtained either in Coulomb gauge or in TAG. Basing 
on this, we will adopt here a pragmatic attitude and use the ansatz Eq. (|32|). It will be 
further explained in section |j. We use imaginary time formalism throughout. 

2 Linear Response Theory 

The potential of a charge Q\ at Xj in the presence of another charge Q2 at x 2 is given 

by 1,| 

V (r) = \llj d 3 * • ^(x) + £ 2 aeff (x) ■ £?(x)) (4) 

where the effective field £± eS created by the non-abelian charge Q1 is the sum of the 
applied field ££ and the induced field S(Ef): 

£r s (x)=m*)+HK) (5) 

Sf is solution of the Gauss' law: 

V-£ a -gf dbc £ b -A c = QZ6 3 (x.-x 1 ), (6) 

where A is the vector potential associated with £. 

Eq. (f|) is manifestly gauge invariant since a trace is being taken over tensor fields 
(or equivalently colour is being summed over) and is therefore a physical potential. 

When we put the charge Q1 in the plasma, the electric field E"(x, t) in the medium 
is coupled to the applied field £° through the coupling Hamiltonian: 

H cxt (t) = J d 3 xEJ(x,t) ■ £?(x) . (7) 

3 Another example of non-Debye screening with a 1/r 10 behaviour in the charge-charge correlation 
has been found in Q 
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From linear response theory, the induced field is 

6(E?{x,t)) = if dt'([H ext (t'),E?(x,t)]) 

J — oo 



dt f J d^£^)([E^,t'),Ef(x,t)]) 



In TAG, 

E«(x,t) = -d Q A«{x,t), (9) 

the correlator of electric fields ([E, E]) can be easily related to the retarded gluon prop- 
agator and the Gauss' law satisfied by the applied field in the presence of a static charge 
Qi, is simply f] 

v-^ = g^ 3 (x- Xl ) , (io) 

where the non-abelian term of Eq. @ vanishes for a static 8 field. The solution is 

tf(x) = -iQl J (^3^ (X " X1) ^ • (11) 

In gauges other than TAG, Gauss' law itself is gauge parameter dependent in a 
parametrized gauge. The non-abelian part of Eq. (|6|) does not vanish, and thus the 
solution S a {x) is not of the same simple form as above and is not necessarily gauge 
parameter independent. In parametrized Coulomb gauge, although Gauss' law is not 
gauge parameter dependent due to A independence of the gauge fixing term, the solu- 
tion is not of the simple form of Eq. (11) because the non-abelian term does not drop 
out. Therefore V(r) is not simply related only to ([El, El]) as in TAG. So that the 
gauge parameter dependence of ([El, El]) shown in |TT], ^Cj] does not mean that V(r) 
is gauge dependent. On the contrary, from the definition in Eq. (|j), V(r) is gauge 
independent. 

In TAG, as a consequence of Eq. (|]), (|[), flTT|) , V(r) is directly related to the 
longitudinal propagator so we have 

V{r) =Q V Q 2 / , , 3 2 . — r (12) 

J (2?r) 3 q 2 + n oo (0,q) 

with r = xi — x 2 . We stress that in other gauges, this advantage of TAG is gone and 
V(r) is related to more than just the two-point function []. 

4 This is the form of Gauss' law in any gauge only at leading order. 

5 It has been argued that in the presence of only one static colour charge, the response of the 
vector potential (A) must point only in that direction in colour space, so that a non-abelian term like 
(A) A (A) vanishes. For this reason, the response field (E) and Gauss' law both become abelian. Using 
this line of argument Rebhan ||, [| argued that the potential of a static colour charge is given by a form 
similar to Eq. (|lj) . This is actually flawed because the response field (E) cannot be expressed directly 
in terms of the response vector field (A) . The response field is actually (E) = (OA) — (OA) + (A A A), 
i.e. the non-abelian term should be the average of the cross product and not the cross product of the 
average, so it does not vanish. Likewise, the non-abelian term in Gauss' law cannot be zero. 
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After going through the usual steps, we get for the colour singlet quark- ant iquark 
potential 0, 

VI \- — ~ 1 9 2 i T ) [°° ^gsingr 
[T) " 2iV 2vr 2 r io g 2 + n 00 (0,g) ' 1 J 

From now on, Iloo is the temperature dependent part of Hoo- 

If Il 00 (0,g) is not even in q as found in J7], [L3|], we have to evaluate the integral 
by closing the contour in the first quadrant of the complex q plane. Contributions 
from possible poles or branch cuts in the first quadrant away from the real axis are 
exponentially suppressed as r — > oo. Eq. (|13D can be rewritten as 

u N 2 -lg 2 (T) T r°° zdze- rz 
^(r>m _1 )^ — / — — - — . (14) 



2N 2n 2 r Jo -z 2 + n 00 (0, q = iz 
Eq. flT4] ) shows that we only need the small q limit of n o(0, q). 



3 Braaten and Pisarski Resummation in Temporal 
Axial Gauge 

In TAG, the gluon propagator is 

D 00 (p) = , D 0l (p) = 

n ( \ - 1 ( r _ PiPi] _ 1 P 2 PiPj h 

lAV) ~ v 2 -G{v)Y 3 P 2 J P 2 -F(p) Po 2 p 2 1 &J 

from which the effective or resummed propagator *D llu (p) can be constructed by replac- 
ing G(p) and F(p) with the corresponding parts of the hard thermal loop (HTL) of the 
gluon self energy, the latter being given by |L4] 



<m£(p) = -8W(27r)- 3 (^ - 4n5 M^p- (16) 

where 



^=4^-— (17) 

and Q = (1,Q). From this, with the Debye mass m 2 = g 2 T 2 N/3, we have for the 
effective *G and *F quantities 

OToo = *F(JP) = " J (P)) , ( 18 ) 
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with 



and 



(19) 



2/2 2 

m / p p 



tIf-? /w J- (20) 



Then *Djj is 



U ^V) p 2_* G ( p )l^ p 2 J p 2 + m 2 (l-/(p)) po 2 ' 



To calculate the gluon self-energy, with two Feynman diagrams, we need both the 3- 
gluon and the 4-gluon effective vertices. From f^, 



poQiQj 
47r p-Q(p + q)-Q, 

Here we have set the external energy, go ; to zero for the static self-energy. Because of 
Eq. (fHf) , this is the only required 3-gluon effective vertex component. For the 4-gluon 
effective vertex, the six permutations of the gluon legs of the 4-gluon hard thermal loop 



can be reduced to just a sum of two terms using identities given in [|14[]. For static self- 
energy and for momentum configuration of the tadpole, these two terms give identical 
contributions so for our purpose, the effective 4-gluon vertex is 



rtM -p, g, -g) = (w) 2 f adc f bcd -2 9ij + 4m 2 / ^ - ™ \ a j • (23) 



dVt PoQiQj 



4 Leading Momentum Dependence of Iloo(0, q) 

The leading term in the static self-energy is well known and is the Debye mass squared. 
To find the leading correction and the momentum dependent part, it is sufficient to 
consider only the zero mode in the internal loop energy [^] of the self energy. This is 
easy to see by power counting of soft momenta. Take for example the gluon self-energy 
graph with two 3-gluon vertices. The usual thermal propagator in Euclidean space 



(27mT) 2 + p 2 



6 



means that for soft momentum and in the zero energy mode it goes as 1/ (gT) 2 while for 
non-zero modes it is 1/T 2 . For a soft loop self-energy graph, we then have l/(gT) 4 for 
two zero mode propagators, {gT) 2 for momenta from the two vertices, (gT) 4 from the 
loop 4- momentum, g 2 from the vertices and 1/g from the distribution function. This 
yields g(gT) 2 , which is g times the leading term. Similar soft power counting can be 
done to establish that there is a contribution of the same order from the tadpole graph. 
For hard loop, there appears to be g(gT) 2 terms left over from the leading hard thermal 
loop (gT) 2 terms but they are proportional to q-Q/|Q|, where Q is the hard internal 
loop momentum vector, and so vanish by virtue of the angular integration. 

With the effective quantities of the previous section, the leading correction is 



AIlg(0,q) = {AllS>,q) + AIlS>,q)}| . , (25) 

l. J zero mode 



with 

(26) 

Anr fe (o,q) = \tY. I ■£^T& d {p,-P,q,-qrD ij (p) . (27) 

As mentioned in the end of section we need only the small q limit of the static 
self-energy when looking for a possible odd term in |q| ]7|, |I3| . This simplifies our 
problem considerably. In fact we can drop the bare vertex part of Eq. (|23| ) because its 
contribution is independent of q. 

We note that since, on the one hand, the effective vertex Eq. (22 ) and the hard 



thermal loop part of Eq. (|23|) are proportional to po and, on the other hand, the 
transverse part of the gluon propagator has no singular factor in po, the zero mode part 
of the transverse-transverse gluon contribution of Eq. fl2~6|) vanishes. The same is true 
for the transverse gluon contribution of Eq. fl27D . 

For the transverse-longitudinal gluon contribution, it turns out that, in the zero 
mode, the bare-HTL vertex and the HTL-HTL vertex combinations are also zero because 



and 




0, (2* 



, (29) 



1 

4vr 4tt Vp-QiP-Q 2 P 2 

where we integrate over the directions of Q« through dfli, i = 1,2. The remaining 
bare-bare vertex combination does not give any odd power of |q|: the integral in p 
reads 

d3 P ( I \2 2, n2\ 1 



(2tt) 3 v " v " ' l ' * yi ' " ' (p 2 ) 2 ((p + q) z + m 2 ) 



P-(P + q) Z - P 2 (P + q) , ,„„ , , 2 , » ■ (30) 
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Because the denominator of Eq. (|30|) is protected by the mass from infra-red divergence, 
we can expand at small |q| < m. It is clear that odd powers of |q| always come with 
odd powers of cos 8, where 9 is the angle between q and p, therefore no odd power of 
|q| appears in the small |q| limit of the longitudinal-transverse contribution. 

Before discussing the longitudinal-longitudinal contribution, it must be mentioned 
that the Matsubara frequency sum is now ill defined for certain terms because of the 
presence of two longitudinal propagators: each has a 1/p 2 , factor and we have only a 
factor of po from each vertex. Therefore terms which do not have any extra factor of 
Po, i > 2 in the numerator are singular. 

In these cases, as was done in || || , we will perform the frequency sums by contour 
integrals using the usual equation [|TJ: 

1 pioQ pioc~\~€. 1 

T^fM = -{ljzf {z)+ l^( m+ f (-;))— J (3!) 

This equation is valid only if f(po) has no singularities on the imaginary axis. This 
is not the case if f(po) contains l/pl or 1/po factors. We therefore employ the usual 
ansatz, namely, that the sum is still given by the above integral. Thus the 1/po and 
1/Pq poles, if existing on their own with no other poles, are vanishing. 

T]T l/p = , T£l/p§ = 0;. (32) 

Let us point out that if one takes the finite T axial gauge propagator directly from 
zero T (i.e. that is one ignores the problem with periodicity of the longitudinal propa- 
gator at finite T |TJ|, see comment at the end of section |5[) then this is the only known 
way at finite T to treat the 1 jp\ singularities. The alternative suggested by Leibbrandt 
and Staley [IT] amounts to performing the uniform replacement 

1 v Po . n 
— — > lim , e > (J 

Po Po + %e 

1 v 2 

-> lim ^ — ? , e > (33) 

Po ^{pl + ief 

and dropping the zero mode. 

The only justification of this prescription is that it leads to the correct answer for the 
leading order Debye mass. This is however not surprising since this calculation is only 
sensitive to hard loop momenta. Furthermore, it is known that due to the plasmon effect 
Jl|], the leading correction to m 2 starts at 0(g 3 ) and not at 0(g A ) as is expected in naive 
perturbation theory. We know by using soft power counting that the 0(g 3 ) correction 
comes only from the zero mode, therefore we argue that any prescription which removes 
the zero mode cannot be the correct one. 
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We calculate the longitudinal- longitudinal contribution by expanding the longitudi- 
nal propagator, making use of m 2 / (p 2 +m 2 ) < 1. We obtain from Eq. ( |21| ) the expanded 
form 

The longitudinal-longitudinal contribution is then 

d 3 p 1 



Ang? L (0,q) = 2p 2 iVT(J a6 



(27r) 3 (p 2 + m 2 )((p + q) 2 + m 2 
1 ( m 2 I{p) \ ni ( m 2 I(p + q) 

ni=0n2=0 PO 



y y y — 

z-^ Z-/ ^2 I p2 _|_ m 2 J Up -|- q)2 _|_ m 2 



|p-(p + q) + m 2 - m 2 (I(p) + J(p + g)) + m 2 p 2 Q J(p, p + g)| (35) 
where J(p) is given in Eq. ( p~9| ) and J(p, q) is 

JM = /l^^WT (36) 

(for more details of this function, see [II]). As the function /(p) is proportional to po 5 
all frequency sums in Eq. ( [35|) with ni +^2 > 2 have vanishing zero mode contributions. 
We are left with the following terms in Eq. fl3"5"|): 

• one frequency sum with a singular zero mode proportional to 1/pq, which vanishes 
according to the ansatz Eq. fl3~2~|). 



frequency sums with zero mode proportional to 1/po- These are performed using 
Eq. (|3l|) and Eq. (|3~2"D. The Bose-Einstein distribution function is approximated 
N(p) ~ T/p for p << T, in order to extract the leading contributions. These, for 
the same reason as Eq. ([30]), yield no odd power of q in the p-integrals. 

frequency sums with well-defined zero modes: we perform the angular integrals 
first and then take the zero mode. The only odd terms in |q| originate from the 
product I(p)I(p + q) proportional to 

go+jpj \ i„ ( go+lp+qj 



dflt dQ 2 1 1 ) lypo-lp+ql ) — > 7T 2 1 

4vr 4tt p .Q 1 ( p + q ).Q 2 4|p||p + q| po=0 4|p||p + q|' 

Collecting the latters and adding an analogous contribution coming from the longitudinal 
part of the tadpole diagram Eq. (|27|), we get 

g 2 NTm 4 ir 2 f d 3 p 1 



3 |p||p + q|((p + q) 2 + m 2 ) 2 (p 2 + m 2 ) 2 

xf^ + q 2 (p 2 + ^ 2 )-(q-p) 2 l • (3* 
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From this we find no linear term but a cubic term in |q| which is g 2 NT\q\ 3 /8m 2 . So 
that, in the small |q| limit, we get the following behaviour for Iloo: 



r tt m \ 2 . 2/tii i i o9 TN( l , 9 TN\q\ , 
hm n 00 (0,g) = m + a£fT|q|+/? h7 5 h 0(|q| ) (39) 



|q|-»0 " " " * m III 2 



where a, (3 and 7 are numerical factors with a = 0,7 = l/8. The factor f3 is not 
as important in front of the q 2 term already present in the potential Eq. (|T3|) since 
it is smaller by a factor of g. Had the factor Eq. (|37| ) not been present, no odd |q| 
term would have appeared. The positive cubic term in Eq. (|39"D is responsible for the 



algebraic decrease of the potential V(r) of Eq. (lj): 



V(r _ x) H (40) 

which behaves as ~ 1/r 6 in agreement with J7|. The presence of this cubic term means 
that the screening potential is no longer dominated by a pole but by some cut starting 
at |q| = 0. On a physical level, the self-energy in TAG is very different from that in 
covariant gauge: we get a longitudinal-longitudinal contribution whereas in covariant 
gauge, this contribution is zero. We find that we have no need to include the magnetic 
mass since transverse gluon does not yield any odd term in the zero mode. 



5 Discussion and Conclusion 

In this paper, using Braaten and Pisarski resummed perturbation theory we study the 
leading correction to the Debye mass in TAG. Indeed, it is the only gauge where the 
potential between two charges is unambiguously related to the two-point function. We 
find that provided there is no branch cut touching the real axis in the first quadrant 
of the complex q plane, one only requires the time-time component of the gluon self- 
energy at small q. This gives a cubic term which yields a repulsive power screening 
potential for the non-abelian theory. This is in direct contrast with the accepted notion 
of Debye mass. Nevertheless, one has to be careful in trusting such a result. We have 
merely looked at the small |q| behaviour of IIoo(0, q) for the dominant behaviour of the 
potential. We expect the same pole as found by Rebhan []3], f| to be present in this 
gauge due to the theorem |18[] that the poles of the propagator are gauge invariant. But 



in this case, the behaviour of the potential is not dominated by this pole and therefore 
it could not be interpreted as the Debye mass. 

Attempts to understand the difference of our result with other gauges may lead one 
to question the use of the ansatz summarised in Eq. (^3) which is known to be true 
up to subleading order in the non-resummed theory || and is untested at higher order. 
However, notice that the crucial cubic term in Hoo comes from frequency sum with 
well-defined zero mode terms. 
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Also, James and Landshoff |l6j have pointed out that the TAG imaginary time longi- 
tudinal propagator should not be periodic in time as in the usual Matsubara framework 
used in this paper: the requirement to include only physical states in the thermal a- 
verage breaks periodicity. They have initiated a rather complex formalism appropriate 
to this gauge based on this. Their new longitudinal propagator has the advantage that 
it is free of po = pole and therefore one does not have to use any pragmatic ansatz. 
Calculating the correction to the Debye mass in this formalism is necessary to confirm 
our result. 

After finishing this work, we noticed Braaten and Nieto |19] and Rebhan pt| have 
extracted the Debye mass from the gauge invariant correlator of two Polyakov loops. 
They found the expected exponential behaviour of the potential and no power screening. 
We do not know exactly where is the source of this disagreement. It could be due to 
one of the problems that we have already discussed above. We leave this as further 
investigation. 
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